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Being motivated by recent experimental data on YBaCuO, we calculate dc resistivity p in the spin- 
gap state of charge-spin-separated t-J model by using a massive gauge theory of holons and spinous. 
The result shows p(T) deviates downward from the T-linear behavior below the spin-gap on-set 
temperature T* as p{T) cx T{1 - c(r* - T^} where the mean field value of d is 1/2. To achieve 
smooth deviation from the T-linear behavior, one needs d > 1. The deviation becomes reduced with 
increasing hole doping. 
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Many quasi-two-dimensional cuprates with high super- 
conducting transition temperature Tc exhibit anomalous 
metallic behavior above Tc in Hall coefhcient, magnetic 
susceptibility, etc. which calls for a new theoretical 
explanation, probably in a framework beyond the con- 
ventional Fermi-liquid theory. Anderson |^ pointed out 
that charge-spin separation (CSS) phenomenon may be 
a possible explanation. 

In the t-J model of strongly-correlated electrons, the 
CSS is naturally described by the slave-boson (SB) (or 
slave-fermion) mean field theory (MFT). When one incor- 
porates fluctuations of MP's, they behave as gauge fields 
coupled to holons and spinous and the system possesses 
a U(l) gauge symmetry. The CSS can be interpreted as 
a deconfinement phenomenon of this gauge theory. The 
system is expected to have a confinement-deconfinement 
(CD) phase transition and the transition temperature 
"^CSS can be estimated 1^, below which the CSS takes 
place. 

The observed T-linear behavior of dc resistivity p{T) 
has been taken as one of signals of the universal 
anomaly in the metallic phase of high- Tc cuprates. Lee 
and Nagaosa Q showed p{T) cx T for fermions and 
bosons interacting with a massless gauge field. This sys- 
tem has some relation to the uniform RVB MFT of the 
t-J model in the SB representation. 

Recent experiments on YBaCuO by Ito et al. ||] and 
others reported that p{T) deviates downward from the 
T-linear behavior below certain temperature T*(> Tc). 
This T* coincides with the temperature determined by 
NMR and neutron experiments at which a spin gap 
starts to develop. So it is quite interesting to calculate p 
in the spin-gap state of the SB t-J model. 

The eff'ective theory used in Ref. |^ is inadequate for 
this purpose, since it assumes no spin gap, containing 
only one gauge field associated with the hopping ampli- 



tudes of holons and spinous. In Ref. 0, we introduced 
multiple gauge fields, and argued that, when a spin gap 
develops, these gauge fields become massive due to a 
gauge version of Anderson-Higgs mechanism. Below we 
use an effective gauge field theory of holons and spinous 
emerged from these considerations and calculate p for the 
region Tc < T < T* to obtain p cx T{l-c(T*-T)'*}, near 
T* where d is the critical exponent of gauge-boson mass, 
rriA oc {T*—TY. This reduction reflects the fact that m\ 
suppresses gauge-field fluctuations and inelastic scatter- 
ings between gauge bosons and holons and spinous. The 
MFT gives d = 1/2, but to obtain a more reliable p, one 
needs to calculate d by renormalization-group method. 
The SB t-J Hamiltonian is given by 
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fxa- is the fermionic spinon operator with spin o'(=t,J,) 
at site a; of a 2d lattice and is the bosonic holon 
operator. The direction index i{— 1,2) is used also as 
unit vectors. 

The MFT is obtained by decoupling both t and J 
terms. For dcfinitcness we follow Ubbens and Lee 

Hmv^Y.{y^Xx.\' + ^\XxA 

- \ E - /li/I+n) + H-c.} . (2) 

The partition function Z(/3) [/? = (fcsT)"!] in path- 
integral formalism is given by integrating out two MF's, 
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the hopping ampUtude, Xxi, and RVB amphtude, \xi de- 
fined on the hnk [x, x + i) |TQ]; 
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I dT[-J2iblk + J2flU)-HMF}. (3) 

In the CSS state (T < Tcss), iXxz) + 0. The spin- 
gap state may reahzes in CSS and is characterized by 
a condensation of A^^^, {\xi) 7^ 0. Let us parameterize 
Xxi = xUxu Uxt = ex.p{iAxi), Xxi = (-)*AT4i, V^i = 
exp{iBxi), assuming uniform RVB. If one ignores phase 
fluctuations by setting A^i = B^i = 0, spinon excita- 
tions has the energy E'^{k) = {(3Jx/4) ^-cosfcj- 
fip}^ + {A J2ii~y h}'^, where /i^? is the chemical po- 
tential to ensure {fl^fx^ + /Ix/xi) = 1 - ^ {S is doping). 
We introduce also for {bl.bx) — <5. There appears a 
spin gap A(cosfci — cosfc2). 

Under the gauge transformation, 6^ = eKp(i9x)bx, 
f'xg. — exp{i6x)fxa, the phases of MP's transform as 
A'^i = Axi + 0x- Ox+u B'xi = Bxi + ex + Qx+;- Their 
behavior can be studied by the effective lattice gauge 
theory ylLGT(C^, V') that is obtained by integrating over 
bx and jxcn e.g., by hopping expansion. For T* < T < 
Tcss I ^LGT contains only ?7's interacting via the con- 
ventional gauge couplings like x''t^2:2 C^a;+2i C^Ii- 
Their quadratic terms for Axi show that Axi behaves as 
a massless gauge field. For Tc < T < T* , A develops and 
new couplings like X^A^fJ+j 2 ^^1+2 1^2; 2 C^x 1 are gener- 
ated. This gives rise to a mass term of Axi, 'm\A^^^ with 
m\ = x^A^. Also, the quadratic Bxi are also massive. 

The holon part A^^ of the effective low-energy contin- 
uum field theory Acq is obtained as 
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where Di = di — igAi (we introduced the gauge coupling 



constant g — 1 for convenience), and (27715) 
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(a is the lattice spacing). The spinon part A^g may be 
written as 
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drd^k {AsG(fc)/T(fc, r)fi{-k, r) + H.c.} , (5) 



with (2mi?) ^ = 3Jxa^/8. Here we introduced the con- 
tinuum version of the spin gap 
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where kF{= \/2mF^iF) is the Fermi momentum of 
spinous and kp — ^27r(l — S)/a at T <ti Tp (Tp is the 
Fermi temperature). In (^, we use the renormalized spin 



gap A(T) defined as (A^,;) = (-)U(r), instead of its MF 
value A to take into account the efi'ect of phase fluctua- 
tions of Xxi effectively. Below we calculate the resistivity 
of the system A^g = ^fg + Afff (ll)- 

The propagator of the gauge held, Dij{x,T) = 
{Ai {x, T)Aj (0, 0)) , is generated by fluctuations of spinous 
and holons, i.e., Dij = {Hp + IIs)"-^, where 

'np^B tj{x,T), = ~{Jp,B i{x,T)Jp^B j(0,0))ipi 

+StjS{x)S{T)np^B, (7) 

representing one-particle-irreducible (IPI) diagrams of 
spinon and holon loops. Jpi = {2'mp)^^ ^^{if^difa + 
H.c.} and JBi = {2mB)^^ {ibdib + I1.C.} arc currents cou- 
pled to Ai, and np = (1 — 5)/a? and ub = S/a^. In 
the Coulomb gauge, the propagator at momentum q and 
Matsubara frequency e; = 2-kI/ [3 is written as 



D{q,ei) = {I{p{q,ei) +IiB{q,ei)}-^ 



(8) 



Since we shall need D later in calculating p{T), we 
obtain D below in the random-phase approximation as 
D ~ (H]j -I- H]?)~^. When the spin gap is sufficiently 
small, its effect to H^ is evaluated by perturbation giving 
rise to a mass term as discussed above; 
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We used the relation, /i_F ~ Tmp/mp and vp = kp/mp. 
The superfluid density of spinous is calculated for small 
AsG(fc)/(fcsT) as 
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tan(/) and |A;| = kp. nj| is given by 
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where /s(e) = {exp(/3e) — 1}^^. Eqs.(|9|,|rT|) above are 
obtained for small q{<^ Ti/a). For large q's, they should 
be replaced by anisotropic expressions due to AsG(fe)- 
These anisotropy can be ignored as long as the spin gap 
is sufficiently small. 

From the linear-response theory and loffe-Larkin for- 
mula p^ , the dc conductivity cr(= an — (T22) is ex- 
pressed as 
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aij = lim lim — -Ily (q, -ie), 

e^Q q^O ~ie 
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where 11, TiF,B are response functions of electron, spinon, 
and holon, respectively. So one has — a~j^^ + Up"^ . 

In the spin-gap state, the spinon conductivity diverges 
CT_F ^ oo due to a superflow generated by RVB condensa- 
tion {\xi) ^ 0. This is an analog of the well-known fact 
in the BCS theory that the electron resistivity vanishes 
below Tc due to a superflow generated by Cooper-pair 
condensation. Actually, has the same structure as 
the BCS model. Thus the total resistivity p = cr"^ in 
the spin-gap state is equal to the resistivity of holons, 
p = <J~^ . Effects of spinous to p certainly exist and show 
up through the dressed propagator -D(<7, £;) in calculating 

Now we calculate the response function ITs. By solving 
the Schwinger-Dyson equation approximately following 
the steps in Ref. jl^, we arrive at 
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ASs(9, En! representing diagrams containing self- 
energy of holons, Eb((7, e„), is necessary to keep gauge 
invariance, 

Aj:B{q, £«; <^i) = ^B^iq, <^i){'^B{q, <^n)G%{q, e„) 

-I]B(g,e„ + eOG|(9, £„ + £;)}, (14) 

However, in the perturbative calculation, this combina- 
tion vanishes in the dc limit by the symmetry under sum- 
mations. We expect this term dose not contribute to the 
dc resistivity, and neglect it hereafter. 

^ B{q, £;), representing vertex diagrams, contributes 

to (Ts, 

g X (g' - g) ] ^ g ■ (g' - g) 
|g' - g| J q^ 

xD{\q' - q|,e„, - e„)i?B(q', e„/; e,), (15) 

where q x q' = q^q'y — qyq'^- We set q of Tb in Ub to 
a fixed vector of typical length qB, q% = Sttub [@|- We 
also fix the length of q' of D in (nsl) be g^. 
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We consider the underdoped region, <C n^? ((5 <C 
1), and temperatures around ^ ub/itib- As- 

suming that D(q,ei) in Tb dominates in the region 
near the static limit e; = 0, we use the upper ex- 
pressions in (p|, pl|). In the denominator of D, the dis- 



sipation term, ^Jn/ {2Tr)\ei\/q, y/E 
larger than the term, (7^/(127rr;ty, m, = 
fB{\f^B\)/{2mB), as long as e; ^ 0, since their ratio is 
smah, {q%/m)/{y/n/{qB(3)} 0{y/n^Jri^). 

So the n'-sum is dominated at e„' = e„. Then we get 
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— {sgn(e„ + ei)~ sgn(e„)} + 0(e?) 
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To calculate nBij(0, e;) we insert ( p^ into ( |l3|) and 
do the g-integral and n-sum as in (|l5|) to get 
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where lim(:j_,o C(T) is finite. 

After analytic continuation e; > 
(|To|), we finally obtain the resistivity. 
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For T* < T < Tcss, nf,{T) = and this result repro- 
duces the T- linear behavior of Ref. 1^ . For T near and 
below T*, one expects the behavior |A(r)| cx (T* - T)"*, 
with a critical exponent d, and we have p(T) oc T{1 — 
c(r* — T)"^}. Note also that the downward deviation of 
p{T) from the T-linear behavior is reduced with increas- 
ing the doping 6. In Fig.l, we plot p{T) of ([20|) with 
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various values of d. The MFT value d = 1 /2 is not con- 
sistent with the experiment. To achieve smooth deviation 
from the T- linear behavior, one needs d > 1. This sug- 
gests that fluctuation effect of phases of Xxi is important 
to obtain a realistic curve of p{T). One could produce a 
reliable curve of p(T) by inserting experimental data of 
A(T) into (|2^), but the available experimental data are 
not enough for this purpose. 

The data |^ show that one may fit p in a form Co + 
CiT for T* < T. This implies spinon contribution to p, 
calculated as ap^ oc T'^^'^ /np Q [Q, is negligibly small 
compared with cr^^ cx T/ub due to higher power in T 
and a small coefficient, cr^^ = for T < T* as explained, 
but the discontinuity at T = T* in cr^^ is not observable 
due to its smallness. 

The constant part Co, surviving below T* , may be at- 
tributed to scatterings of charged holons with impurities. 
They may be described by i?imp = J2 ^xbl-bx, where 14 
is a random potential. Actually, standard calculations 
show that it generates T-independent contribution to p, 
AtT^^ cx l/ns at intermediate T's |l5). 

As T goes very low, 11^ is dominated by the mass term, 
Ilp{q,ei)/g'^ — np{T)/mF, while Ilf does not change. 
The Landau damping term from holons is smaller than 
the mass term above. This case has been studied in Ref. 
psf , giving the result cr^^ oc . Thus, in p{T), weak- 
localization effect by impurities, pwL ~ CwLlog(T), 
dominates over inelastic scatterings by gauge bosons. 
This situation is in contrast with the effective gauge 
theory of two-dimensional electrons at half-filled Lan- 
dau level, in which the transverse mode of Chern-Simons 
gauge field remains massless down to T = and renor- 
malizes Cwl 
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FIG. 1. Plot of the resistivity p{T) divided by 
po = 2TvkBT*mF/{e^nF) for A(r) ~ Ao(l - T/T*)'* 
{d = 1/2,1,2). For definiteness we chose 5 — 0.05, 
Ao = 2fcflT*/7r and 3m/(4mF) = 0.5. 
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